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ON POSSIBLE SIMILARITY SOLUTIONS FOR THREE-DIMENSIONAL 
INCQMEKESS IBLE LAMINAR BOUNDARY-LAYER FLOWS OVER 
DEVELOPABLE SURFACES AND WITH PROPORTIONAL 
MAINSTREAM VELOCITY COMPONENTS 1 
By Arthur G. Hansen 


SUMMARY 

This report presents an analysis of possible similarity solutions of 
the three-dimensional, laminar, incompressible, boundary- layer equations 
referred to orthogonal, curvilinear coordinate systems. 

Requirements for the existence of similarity solutions are obtained 
for the following two cases: flow over developable surfaces; flow over 

nondevelopable surfaces with proportional mainstream velocity components. 
The analysis obtains permissible forms of mainstream, velocity components, 
the square of differential of arc length on the surface, and the simi- 
larity parameter. A basic class of surfaces is found from which all other 
permissible surfaces may be obtained. 

Necessary and sufficient conditions are found for expressing the 
ordinary differential equations resulting from the similarity transforma- 
tion in uncoupled form. The analyses shows that uncoupling is possible 
only when the surface is developable and a surface coordinate system 
characterized by (ds)^ = (dx^)^ + (dx^)^ is employed. 


INTRODUCTION 

Theoretical research on the important problem of three-dimensional 
boundary- layer flow has been greatly restricted because of the complex 

-kuhe information herein was originally presented as part of a thesis 
entitled "Similarity Solutions of the Laminar Incompressible Three- 
Dimensional Boundary- Layer Equations" submitted in partial fulfillment 
of the requirements for the degree of Doctor of Philosophy at the Case 
Institute of Technology, May 1958. 
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nature of the governing equations. However , a degree of success lias been 
achieved in incompressible-flow analysis by searching for exact solutions 
of the equations for special types of mainstream flows. These exact so- 
lutions of the three-dimensional, laminar, incompressible boundary- layer 
equations have generally been formulated by use of "similarity” tech- 
niques (e.g. , refs. 1 to 7). In instances where the technique is appli- 
cable, the partial differential' equations of the boundary layer reduce 
to a system of ordinary differential equations. The corresponding solu- 
tions for the boundary- layer velocity components are such that the veloc- 
ity profiles differ by, at most, scale factors along coordinate directions. 

To date, similarity analyses have been carried out for only very 
special types of mainstream flows and for special types of coordinate 
systems. A question has always existed, therefore, as to the most gen- 
eral class of problems amenable to a similarity analysis and as to what 
limitations are inherent in its use. This question Is successfully 
answered for two-dimensional laminar flows in references 8 to 12. These 
investigations show that, in general, similarity solutions can only be 
found for mainstream flows in which the velocity varies as a power of the 
distance along a surface or as an exponential. Recently, several inves- 
tigations have attempted to determine conditions under which similarity 
solutions exist for laminar, incompressible, three-dimensional boundary- 
layer flows (refs. 13 to 17). In investigations of this kind, the 
type of coordinate system employed plays ah important role because of 
similarity of velocity profiles in coordinate directions. Consequently, 
reference 13 considers the special case of a stationary rectangular co- 
ordinate system and determines what possible mainstream flows referred 
to such a system lead to similarity solutions. Reference 14 determines 
permissible mainstream flows referred to polar coordinates; reference 
15 determines mainstream flows confined to regions of small angle vari- 
ation with respect to polar coordinate systems. 

The research presented in reference 16 is more general in scope than 
that in reference 13 or 14. In reference 16 the" coordinate system 
assumed is an arbitrary orthogonal, curvilinear system. In this respect, 
the approach employed therein is similar to that which will be used in 
the present report. However, the complexity of the problem has required 
the application of additional assumptions, both In reference 16 and in 
the analysis given here. In reference 16, for example, the mainstream 
flow is assumed to be irrotational. In the present analysis, two inde- 
pendent assumptions (each different from that of ref. 16) are made and 
investigated separately. 

The .first assumption which Is applied herein Is that the surface 
over which the flow takes place is developable. Examples of such sur- 
faces are cones, cylinders, and, of course,_the plane. Such surfaces 
are characterized by the geometric property of zero Gaussian curvature. 

The second assumption pertains to a restriction on the form of the 
boundary- layer velocity components when the flow takes place over a _ 
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surface which may be nondevelopable . It ■will be shown that the assump- 
tion is equivalent to specifying proportionality of mainstream velocity 
components in reference to a particular coordinate system embedded in 
the flow surface. 

The results of the present study are compared in detail with those 
of reference 16 on page 40. However, it might be well to point out here 
that there is a major distinction between the philosophy of reference 16 
and the present treatment. Here, greater emphasis is given to the geo- 
metric aspects of the problem. In reference 16, a great deal more atten- 
tion is given to establishing necessary and sufficient conditions for the 
existence of similarity solutions. 

In the following sections, general requirements for similarity solu- 
tions will first be derived. Then, further developments on the assump- 
tion mentioned previously will be presented, and specific conditions will 
be determined for reducing the boundary- layer equations to ordinary dif- 
ferential equations. 

Acknowledgement is made to Dr. Gustav Kuerti of Case Institute of 
Technology for his interest and advice in the preparation of this work. 


CONDITIONS FOR SIMILARITY SOLUTIONS 


Consider a surface in space in which an arbitrary orthogonal coordi- 
nate system (x^xg) has been embedded. Let y* be a coordinate normal 

to the surface (fig. l) . The boundary- layer equations referred to such 
a system are the following (refs. 18 or 19): 






+ 



+ UjUgk-L 




U-, du-, u 2 du, „ 

x + + "iWl - u i*2 

3 v du R . - 2. 


U„ 5Uo _ 3u, O d 2 i 

+ — h + _L_ _£ + ut U 2 kc - uTkn 

*2 y; 122 11 &yS 

U-, dUo Up dUo 2 

-i ^- 2 . + U,U„k„ - U7k, 


Ju 2 


’ u 2 k l 


Jo 6Uo 

454 + w 2 “% 

1 Su i i du 2 1 Sv 

^ 55 + s; + 


(la) 


(lb) 


(lo) 
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where 

u l> u 2-> v boundary- layer velocity components in thex-j_-, x 2 -, and y*- 

directions, respectively 

y y*/ >/v 

k ,k 0 geodesic curvatures 8, of coordinate lines *= const, and 

1 dt '-=t dt 

1 ^^1 1 ^2 
X1 = const., respectively, kj. - EJhJ 


h;j_,h 2 metric tensor components related to the differential of arc 
length on surface by (ds) 2 = hpCdxq) 2 + h§(dx 2 ) 2 

Ui,U 2 inviscid main-flow velocity components in vicinity of surface 

(All symbols are defined in appendix A.) The boundary conditions are: 
At y = 0, U]_ s ug = v = 0 

lim Uj_ = Uj lim u 2 = U 2 . 

y -*■ • y -*■ • 


Two possible situations might now be noted relative to the main- 
stream velocity components. The first is when neither Uj_ nor U 2 is 

identically zero. The second is when one of' these components is iden- 
tically zero,, and, consequently, one of the coordinate lines corresponds 
to a streamline of the main flow. For these two cases_(when neither Up 
nor U 2 is identically zero) follow the classical approach for obtain- 


ing similarity solutions , 

and assume that 

Uj_ and u 2 are expressible 

as follows: 

u i 

= ^(T)) 

(2) 


u i 

dT] 


u 2 

-Mni 

~ (3) 


U 2 

dt) 

a The geodesic curvature at a 

point P 

of a curve embedded in a 


surface in space is numerically equal to the curvature of the plane 
curve obtained by projecting the surface curve onto the tangent p lan e 
to the surface at P. ~ 
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where F(t|) and G(t]) are as yet undetermined functions of the similarity 
variable T} defined as 


y* 

n = /= g(x x ,x 2 ) = yg(x 1 ,x 2 ) 

yv 


U) 


and where g(xi,x 2 ) is an arbitrary function of x^_ and x 2 . (See ref. 
13 for a detailed explanation on specification of the form of t\. ) 


For the case of one mainstream velocity component identically zero, 
it can be assumed without loss of generality (from symmetry considera- 
tions) that U 2 s 0 and u-j_ and u 2 are expressible as 


^1 _ dF(n) 
Up - dr; 


( 5 ) 


_^2 dG(ti) 

Uf = dn 

where is a function that will be determined in the analysis. 


( 6 ) 


Transformed Equations for U-l and U 2 Not Identically Zero 

From the continuity equation, equation (lc), and the definitions 
for Uq_ and u 2 , it is possible to obtain an expression for the 

boundary- layer velocity component v. Substituting equations (2) and 
(3) into equation (lc) and using equation (4) give 


UqF'kg 


U-lF" 



P' fcU-L 
+ hq <3xi 


+ UgG f k^ 



+ 


Gi au 2 x 1 bv 

*2 ^2 yfi ' 5 * 


0 


(7) 
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Solving equation (7) for 6v/6y and integrating from 0 to y with 
respect to y result in 


v = - 


ZL 

e 


X d ** 2 + 1 § - § vtf ) + 


G (Ugk-L + — - 2^ - dln_g\ + Hi dinjg -p, + - ^2 9l£j ,1 

V 2 1 h 2 h 2 dx 2 ) * h x Sx x V + h 2 <bx 2 110 J 

fU-^Xg) 


( 8 ) 


where 


f(x 1> x 2 ) = v(x 1 ,x 2 ,0) + 





1 Su l _ U 1 5 In s 
hi 5xi hi ^Xi 


G(0) 


( U 2*l + i“ 




Su 2 _Uz 

Sx^"h 2 


6 In 
bxo 


I 


At this point, it would be well to discuss restrictions on the func- 
tions F(t)) and G(rj) resulting from the boundary conditions on Ui, 

u 2 , and v. From equations (2), (3), and (8): 

(l) The boundary conditions «= ug = 0 for y = 0 imply 


( 2 ) 


( 3 ) 


F'(0) = G'(0) =0 


The 

conditions 

lim Ui = Ui 

and 

lim ug = Up ” imply 



y*-« 


y-x° 


lim F 1 

'(n) = 

1 and 

lim G 

*0i) - i 







The 

condition 

O 

II 

!> 

for ; 

Y = 0 ; 

implies 

"r 

VV 

f(x x ,x 2 ) = 

F(0) 

(“1*2 + 

i au i 
hi Sx£ 

_ Ul ^ 3a &) + 

hi dxi J 



0(0) 

(“2*1 + 

_L* U 2 
h 2 3x2 

u 2 6 In g\ 
h 2 dx 2 / 
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Now, it can be shown by a slight 

in reference 13 that there is no loss of generality \ _ 0 . 

Sat F(o) = 0 and G(0) = 0, which in turn implies * “ °* 

This result will be used in the following development. 

Substitution of equations (2), (3), and (8) into equations (la) 
and (lb) yields, respectively, 


h l ° X 1 V^l 


su x u x a m 


hi 55 2b]_ oxi 


+ Unko FF" - g 2 F'" + 


/U 2 B in U X \ _ (1 ^2 _ JL i^-S! + ^ GF" - 

\hj ^ X2 +U 2 k lJ GF ^55 2hg Sx 2 2 i; 

U§ A 3Ut Uo S in U x u| \ . . 

*2<3 ' 2 - + ^ - 57 - + k2 J = 


1 S0 2 „ 




(I ^ ^ - - fe 5 ■ - k : ^ + **) re " - 

<10) 


^ kiF' 2 


_ + 4^0 rq'l nnfl ( 10 ) will be temed the "transformed equations," and 

»«S“srSiI 

X^Sntica^ero. 
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Transformed Equations for One Mainstream Velocity- 
Component Identically Zero 

How consider the case U 2 e 0 and define uj. aid u 2 by aqua-' 
ObtSnefl 6 ?“ boundary-layer yeiocity component v ” 

^thV°SpSedV^. be laen “° Sl *> ^ 1D - 

Substitution of the expressions for Ul , u 2 , and v into eaua- 
tions (la) and (lb), respectively, then gives 


£~^ F ' 2 - SI + U l k 2 - ^ FF" - g2p... + 

($ ii ? P 1 + & \ . (x Su 2 u 2 a in * \ 

V h 2 ^*2 2 1 ) (h 2 " 2hg <hc 2 K + U 2 k lj GF " 


.*2 


u 2“ 1 dU n 

Ut k 2 G ' "h^^ = 0 


( 11 ) 


2h 2 dx 2 


2 3 In g2 \ P 

®- GG - g^G 1 ” + 


i SU 2 p ,2 /l &J? * 

h 2 ^ ^ + U 2*l 

/UldlnU*^ \ _ A 5Ul U^ aja 2 t 

v h 2 ^ X 1 1 7 (h x 5x^ 2h 1 Sxj_ g + U l k 2j l 1 ®" - 


U 


u; 


uf k i F ' " u| k i 


kn = 0 


( 12 ) 


(9) aS^n? 1X i ? ld ^ 12) are of exactl y the same form as equations ‘ 
(l ? ) + exce ^ t for tke absence of several terms vhlch are ^t Sef 
ficients of terms involving P, G, or their derivative^ OMs fact 

in £ r 1 +hp r + a ^ t ^ de ^ erm±ning the Several conditions necessary for reduc- 
ing the transformed equations to ordinary differential equations. 


i 
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It should he stated that the boundary conditions for G' (t)) for 
this case are given by 

G'(0) = 0 
lim G'Ol) = 0 

T) 

The boundary conditions on F'(t)) are the same as for the previous case. 


Necessary Conditions for Obtaining Similarity Solutions 
co 

bo It is obvious that equations (9) and (10) or equations (ll) and 

° (l2) become ordinary differential equations if the coefficients that 

depend on x-j_ and x 2 are proportional. This assumption will now be 

introduced (hereafter called assumption A) . The assumption of propor- 
tionality is a sufficient condition for the solution of the system of 
equations to depend only on t). Whether this is a necessary condition 
will be discussed on page 40. At present, proportionality of the coeffi- 
cient is assumed simply as an additional hypothesis under which to solve 
. the systems. First, examine the implications of this assumption rela- 

tive to equations (9) and (10) . As the function g 2 serves as a coef- 
ficient in both equations, the general requirements are as follows: 


= a (x dU2 - ^2_ 

4 (h.2 3x2 ” 2h2 


/U 2 d In U x \ 

& 3 (hi ”"3x2 + U 2*l 

A U 2 a In Up up \ n auo 

= a 6 fc^ + ^^r- + U 2 ki -ui k2 ) 

^ a in u 2 \ uf 

+ U l k 2 j ~ a 9 k l 

uf 


/U-, a In U 2 1 au ? u? \ 

a iQ (hi S£- * sej + u i k 2 - uj * 1 ) 


W 



10 


MCA 2M 1437 


Now, the following observations are made. First of all,, g cannot 
be identically zero (see eq. (4)). Therefore, if the individual ter&s 




h l h 2. 


U 2 a In U x 


UZ 


+ and k 2 are each proportional to g‘ 

2 .1 « v 


identically zero), the term 


% 

\ h l ^l h 2 


(or 

5 In IX, ** uf \ 

+ ^ h 2 Will be 


proportional to g 2 (or identically zero) . Similarly, if 
a In U g 

\ho 




°i4 t 


au 2 , 1, 

“ 4 ' : k i 


(or identically zero), the term 



are each proportional to g c 


au-- 


Itf 


+ 3 ^ + U l k 2 “ ^*1 


will be proportional to g 2 (or identically zero) . The conditions im- 
posed by equation (13) therefore reduce to the requirement that nine co- 
efficients be proportional. Furthermore, various terms in certain of 
these coefficients can be deleted when these terms are required to be 
individually proportional to g. This leaves the following list of nine 
terms that must be mutually proportional (or identically zero) if simi- 
larity solutions of equations (9) and (lO) are to be obtained: 


^ 1 SU 1 


U l k 2 


g 


^i a In g 2 
“ * §xV 


2hi 


U P a In IX, 

© 4 + U2ki 


u 2 ki 


u. 

t-» y 

2ho Soj 


2 a In g" 


©^*2 


© s? 


U x a In u 2 




+ % 


<D £*i 

u 2 


y 


(14) 
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Proportionality among the preceding terms is equivalent to specifying a 
system of partial differential equations. Solutions of the equations 
determine mainstream velocity components U-j_ and Ug and the components 

hj_ and hg of the metric tensor associated with the orthogonal coordi- 
nate system. As solutions of these equations will result in equations 
(9) and (10) being reduced to ordinary differential equations, a stated 
proportionality between any two terms will be called "an ordinary differ- 
ential equation condition" and will be abbreviated "o.d.e. condition." 

If o.d.e. conditions are set up for equations (ll) and (12) , it 
will follow that these conditions will be exactly the same as those given 
for equations (9) and (10 ) , except for the replacement of U 2 by U*. 

Hence, solutions for various functions in one case will correspond ex- 
actly to those in the other case. In setting up the equations for the 
o.d.e. conditions, TJg is therefore used to denote either Ug or Ug. 

A general analysis would now involve finding solutions of the sys- 
tem of equations derived from set (14) without further restricting 
assumptions on the nature of the unknown quantities. As mentioned in 
the INTRODUCTION, this most general problem has not as yet been analyzed 
because of inherent complexity. Consequently, the problem will first be 
attacked under the assumption that the flow surface is developable. 

This is equivalent to specifying that the Gaussian curvature K of sur- 
face is identically zero. In turn, this assumption leads to an addi- 
tional equation involving h^ and h^. 

Following the analysis for developable surfaces, the problem will 
be solved under the assumption that U-^/Ug is constant with no restric- 
tion on surface geometry. It will be shown later that a consequence of 
this assumption (for Ug. = Ug) is that streamlines of the main-flow cross 

coordinates lines at a constant angle. 

The requirements for similarity solutions will be derived according 
to the following topical scheme. (No attempt will be made here to solve 
the corresponding ordinary differential equations that result.) 

(l) Hypothesis: K = 0 (pp. 12-24) 

(a) k£ + k| / 0 

(b) = 0, kg = 0 


W 
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(2) Hypothesis: = cUg (pp. 25-38) 

(a) kjkg ^ 0 

1. kj_ and kg nonconstant 

2. kj_ and kg constant 

(k^kg ^ 0 and only one choice of k^ constant will 
be shown to be impossible) 

(b) k^ = 0 ;; - 

ANALYSIS OF POSSIBLE SIMILARITY SOLUTIONS FOR' FLOW OVER 

SURFACES ISOMETRIC WITH THE EUCLIDEAN PLANE 

In the following sections it is assumed that the flow takes place 
over a developable surface (i.e., a surface isometric with the Euclidean 
plane). Relative to this assumption, solutions for h-j_, hg, k-j_, kg, U^, 

U 2 , and g^ will be obtained from the o.d.e. conditions. Before this 

can be carried out, however, certain relations are needed from the in- 
trinsic geometry of surfaces . These relations will now be developed. 


Coordinate Curvature Relations for a Surface of 


Zero Gaussian Curvature 


Two surfaces are said to be isometric if it is possible to find a 
coordinate system embedded in one surface which has the same surface 
metric tensor components as a coordinate system embedded in the other 
surface. Furthermore, if two surfaces are isometric, an invariant ex- 
ists called the "total" or "Gaussian" curvature, which will be the same 
for both surfaces. For a system of orthogonal coordinates it can be 
shown that the expression for the Gaussian curvature is 


1 r> 

_ ( 1 


1 5h !\ 


1 \ h l h 2 “1/ 

^2 V 

h l h 2 ^2 /_ 


(See ref. 20, p. 169.) 


(15) 


Substituting the 
equation (15) gives 


expressions for coordinate geodesic curvature into 


K 


■-fe 


dkg 

dxi 


_1_ 

*2 


dk 

dxp 


{)- * 


4 ) 


A. 


*1 


(16) 
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in 

r- 

o> 

■<# 


Now employ a theorem from the geometry of surfaces which states that a 
necessary and sufficient condition that a surface he isometric with the 
Euclidean plane is that the Gaussian curvature he zero (ref. 20, p. 168). 
Hence, for the case under consideration, from equation (16), 


T La , <- 

k, + k<- 


"fe 





(17) 


With this condition relating the coordinate curvatures and metric tensor 
components, and with the o.d.e. conditions referred to in the previous 
section, it is possible to determine the permissible coordinate systems 
and mainstream flows leading to similarity solutions of the houndary- 
layer equations. Equation (l7) clearly shows that k-j_ and kg cannot 
he nonzero constants. 


In general, it will he necessary to make certain initial assump- 
tions regarding the numbered terms appearing in set (14) before- a unique 
set of o.d.e. conditions can he considered for analysis. This follows 
from the fact that an equation expressing proportionality between two 
terms can only make sense if the terms are not identically zero. Con- 
sideration is first given to coordinate systems in which kj_ ^ 0 and 
k 2 ^ 0. 


Reduction of O.D.E, Conditions for 0, k 2 ^ 0 

Under the assumption that k-^ ^ 0 and k^ ^ 0, the terms numbered 
(§ y (?) 9 and (?) in set (14) cannot be identically zero. Hence^ the 
following o.d.e. conditions prevail: 


2 2 

u“ k 2 - a x g f 0 

1 

U l 

% h. - a 2 g a 2 + 0 


(18) 


(19) 


w 
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Differentiating equation 


(18) with respect .to 


X]_ gives 


a + s! . 2 if! si M 

SI i \u§ ug SJ 3xj) 


= a. 


2 ( jin .g 

ui 1 V W 


a IdU, 

+ — g— 


- 2 



2 \-^r * szr 


a In Uo\ 

2 ~srJ 


( 20 ) 


Similarly, by differentiating equation (19) with respect to xg it 
follows that 


Sk l V 


f 


a In g £ 

3£T 


a In Uo 

+ “^T": 


ainui\ 

Sx 2 / 


( 21 ) 


It will now be shown that the right sides of equations (20) and 
(21) can be simplified in form. Following the simplification, equa- 
tions (20) and (21) will be substituted into equation ( 17), and a rela- 
tion between U-j_ and U 2 will be established. _ 


The assumed proportionality between the coefficients listed in set 
(14) leads to the following set of statements: 

(1) From the o.d.e. condition for (2) and @ , 

: 

a in U 1 ' _ 

(2) Either — — s 0 f or from the o.d.e. condition for (l) 

and (D f 


3 In 0 X g 2 h^_ 

“3*1 = a ± ~ 
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(3) From the o.d.e. condition for ® and 


-2 


d In U 2 

Sxi 


* kp + 


a 5& 

Ui 


2hn 


(4) From the o.d.e. condition for (3) a.nri 


a - fe + 




a 6 g 


1 U 1 2h 2 
2 


S laU, 


(5) Either ^ 2 0^ or., as an o.d.e. condition for 

and (7 ) , 


a In U 2 

87 ”7? 


g 2 hg 


Up 


(6) From the o.d.e. condition for (3) and (4 ) , 


a In U, 


-2 


ax. 




From these observations it follows that equations (20) and (21) must 
be expressible as follows: 


at. 


cix£ - 2k x h 2 (nk 2 + a 9 j 

aki / _2) 

^ = 2k 1 h 2 (mk 1 + a 10 ^ 


( 22 ) 

( 23 ) 


where 


n = 0^ 1, or 2 
m = 0, 1, or 2 
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As a consequence of the hypothesis k^_ 0, k g ^ at least one of the 

quantities n, m, a g , or a-j_ 0 must he nonzero. Now/ by using equations 
(22), (23), (18), and (19), equation (17) can be rewritten as 


- 1 . 


k l 1 1 + 2m + a 10 &2 =5 / + k 2 U + 2n + a-|_ a g — I = 0 


U; 


Uf 


(24) 


An o.d.e. condition between @ and is 


u| 




Ul k 2 = wT kl 


'U 2 \ s 

Un‘ 


= a 


ii e; 


(25) 


Substituting equation (25) into equation (24) results in 


( 1 4" 2m 4* a io a 2 ^ 


*11 


— 1 

) 4 - 1 4 - 2n 4 - a^ a^H = 0 


(26) 


where 



From equation (26) and the definition of H, it follows that H must be 
a positive constant. Therefore, the following lemma is established. 

Lemma 1 . - Given an orthogonal coordinate system in which neither 
k}_ nor k£ is identically zero, a necessary condition for equations 

(9) and (10) or equations (ll) and (12) to possess similarity solutions 
under assumption A is that 

Ui = c/G 2 (c-l a constant) 

where Ug = Ug w ^ en considering equations (9) and (10) and Ug = U* 
when considering equations (ll) and (12) • 

Furthermore, from lemma 1 and equation (25) there follows lemma 2. 


r " 




4975 



*±^ IO 


NACA TO 1437 


17 


Lemma 2 . - Given an orthogonal coordinate system in which neither 
kq nor k2 is identically zero, a necessary condition for equations 
(9) and (10) and equations (ll) and (12) to possess similarity solutions 
under assumption A is that 

kl *» Cgkg (eg a constant) 

With the establishment of lemmas 1 and 2, it now becomes possible 
to determine explicit forms for the coordinate - line curvature kq and 

kg the functions h]_ and hg through solutions of o.d.e. conditions 

from set (14) . Hie details of the analysis are presented in appendix B. 
The results obtained can be summarized in the following theorem. 

Theorem 1 . - Given an orthogonal, curvilinear coordinate system em- 
bedded in a surface of zero Gaussian curvature and in which neither k^ 

nor kg is identically zero and neither nor Ug is identically 

w zero. Then, necessary and sufficient conditions for equations (9) and 

V (10) or equations (ll) and (12) to be reduced to ordinary differential 

equations under a similarity transformation and assumption A are that 


(l) h-L = apj^x^qgCxg) 

(a ^ 0) p£( Xl ) # 0 

h 2 = bp 1 (x 1 )q-(x 2 ) 

(b { o) q.g(x 2 ) # 0 

(2) “ 7 V 7 T 

1 bp 1 (x 1 )qg(xg) 


L 1 
2 ap 1 (x 1 )qg(x 2 ) 


(3) g 2 = C 1 p5 _1 (x 1 )q2" 1 (x 2 ) 

^ 0 

(4) = C 2 pJ(x 1 )q“(x 2 ) 

a, t 0 

U 2 = C 3- P l^ x l^ <1 2^ ;x: 2^ 

C 3 / 0 

Equations (9) and (10) ‘becamej respectively ^ 


(5) AjF ' 2 - A 2 FF" - A 3 F ,n + A 4 G'F' - A^-GF" - AgG ' 2 - (A x + A^ - A-} = 0 
B^G' 2 - A 5 GG" - A 3 G ,m + B 2 G*F' - AgFG" - B^’ 2 - (B^^ + Bg + B 3 ) = 0 
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where 



(n + 3)Cg C^n 

A 2 2S B 1 = “ 


C ? (n + 1) 

a 3 - Cj. b 2 = 

C- c§ 7 

A 4 - - On + 1) B 3 = — 

A s = || (m + 3) 


Equations (ll) and (12) reduce to, respectively, the same forms as those 
given in (5), except that the first equation has only Ag as a constant 
term and the second has only B 3 as a constant term. 


Coordinate Systems for Which k^_ ^ 0, k 2 ^ 0 

The following two questions are now posed: 

(1) Is there a unique coordinate system corresponding to the metric 
components h]_ and h 2 defined by condition (l) of Theorem 1 ? 

(2) What is the nature of coordinate systems defined by these 

equations ? _ ...... 

The approach to the previous questions can be simplified somewhat 
by carrying out the following transformation of coordinates. Let 


X x = Pi(x x ) 

x 2 = <3.2 ( x 2^ 


(27) 


Hie square of the differential of arc length in the (X^,X 2 ) system is 
then given by 


(ds) 2 = a^CdX^ 2 + b 2 X 2 (dXg) ‘ 


(28) 
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The coordinate lines X^_ = constant will coincide with X]_ = constant, 
and the lines X 2 = constant will coincide with x 2 = constant. The 
metric components are given by 


h^ = aX 2 

4 = 

The principal difference between the (X-^Xg) system and the (x^,x 2 ) 

system will be in the measurement of length along the coordinate lines. 
A transformation of the previous type will hereafter be called a "change 
of scale" transformation. 



The surface in which our coordinate system is embedded is assumed 
to be isometric with the Euclidean plane. Hence, the question of whether 
or not a unique coordinate system exists with metric components defined 
by equations (29) is equivalent to asking whether or not a unique trans- 
formation of variables exists of the form 


Y 1 = *l(Xi,X 2 ) 

Y 2 = YgCX^Xg) 


(30) 


where and Yg are Cartesian coordinates; that is, 

(ds) 2 = (dY x ) 2 + (dYg) 2 


This question is discussed in reference 20. It is shown in this refer- 
ence that a necessary and sufficient condition for a suitable transfor- 
mation of the form (30) to exist is that the Riemann-Christoffel tensor 
formed from the h* be a zero tensor. For a surface having zero 

Gaussian curvature, this condition is satisfied. Furthermore, the 
transformation will define a unique rectangular coordinate system ex- 
cept for possible translations and rotations. 


The second question on the nature of the coordinate system requires 
a more detailed investigation, which is presented in appendix C. It is 
shown there that the coordinate lines became logarithmic spirals when 
referred to the Euclidean plane as a developable surface. If p and 6 
denote polar coordinates, the equations for the coordinate lines 
Xj_ = X° « constant and X 2 83 X§ *= constant are given, respectively, by 


P 

P 


= c(X°) d e" 

l/d 2 e/d 


■0d 


=C(X£) 


(C8) 

(C9) 
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where 


2 _ a 2 b 2 
° = a 2 +b 2 

2 a 2 
d^ = 


A typical network of such coordinate lines is shown in figure 2. 

Permissible main-flow streamline shapes in spiral coordinate syste m 
(k]_ ^ 0 } kg ’f 0). - The equations for main-flow streamlines in the 
(x^,x 2 ) coordinate system or equivalently in the (X^,X 2 ) system are so- 
lutions of the equation 


Un 


hi dxi 


hT dXi 


aX 2 dX-L 


U 2 h 2 dx 2 h* dX 2 bX l dX 2 


From lemma 1 , U^/O^ = c^. Hence } 


aX 2 dX-L 
bX^ dX^ 


= C 1 


Equation (32) has the solution 

Xj_ = (const.) X 2 


( 31 ) 


(32) 


(33) 


where 


e,b 

r = -A- A o 
a 


This, then, is the defining equation for the permissible main-flow stream- 
lines in the (Xj_,X 2 ) coordinate system. Once again, a more familiar form 

for equation (33) can be obtained by expressing the equation in terms of 
(p,S) coordinates. The transformation between the (p, 0) and (X 1; X 2 ) co- 
ordinate systems is given in appendix C by the equations 


p = cXiX 2 


0 = d In Xt - 


in X 2 


(C5) 
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Substituting equation (33) into (C5) gives 

■y, | 1 

p = PqX7> (Pq a constant) 

0 = ^dr - In Xg + 0 q (0q a constant) 

If Xg is eliminated from equations (34), 


(34) 



e m(0-0 o ) 


(35) 


where 


d(r + 1) 
2 r - 1 


(d 2 r j 1) 


If m ^ 0, equation (35) is the equation for families of logarithmic 
spirals in a planar polar coordinate system. If m = 0, there are 
circles in such a system. Finally, if drr = 1, equation ^34) gives 
9 = constant, that is, radial lines in a planar polar coordinate system. 

Transformation of coordi nates to a basic system . - Before leaving 

the analysis of this system, both the form for (ds)^ and the equation 
of the main-flow stre aml ines will be put in a particular form for later 
reference. First, a scale transformation is introduced: 


X 


1 



X 


2 


x 2 /b 


Equation (28) then becomes 


(as ) 2 




(36) 


The equation for U-|_ = cgUg "becomes 


U 1 - 


n^ + m‘X 2 


e 


(37) 
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In the (X]_,X 2 ) coordinate system the equation of the main-flow stream- 
lines assumes a particularly simple form. As the metric tensor compo- 
nents are identical, the governing differential equation is 

£2 = ^2 

dX x ” 

Hence, the streamline equation is 

a 'X-l + h 'X 2 + c ' » 0 

It is of interest to note that the streamlines in the various co - 
ordinate systems previously employed cross the coordinate lines at con - 
stant angles . This follows from an examination of the formula for the 
cosine of the angle 9 between a streamline and a coordinate line. For 
flows in which the mainstream velocity component U 2 0 and * c^L^, 

(h-|_ dx^)/(h 2 dx 2 ) = c^ on a streamline in a coordinate system (x-j_,x 2 ). 
The cosine of the angle between the streamline and an" x^_ coordinate 
line is then (ref. 20, p. 150) 

2 d*! dx£ “ 2. 

cos 9 - h x dx x ) cIx-l) (1 + c |) ~ 1 + c| 

9 = constant 

A similar result holds for the angle between the streamline and an x 2 - 
coordinate line. 



Seduction of O.D.E. Conditions When One Coordinate 
Curvature is Identically Zero 

Wow solutions of the o.d.e. conditions are considered when the 
curvature of one set of coordinate lines vanishes while the other cur- 
vature does not identically vanish. Initially, it is assumed that 
hi = 0 and k 2 0. The case k 2 = 0, kj_ ^ 0 will follow directly 

from symmetry considerations. 

The analysis based on the o.d.e. conditions and K = 0 leads 
directly .to the following theorem, which is proved in appendix D for 
kj_ = 0, k 2 0. — 


Theorem 2 . - Let (x]_,x 2 ) be a curvilinear orthogonal coordinate 
system embedded in a surface of zero Gaussian curvature and in which 
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in 

r- 

<j> 


k 2 is zero and is not identically zero. Then, a necessary condi- 

tion for equations (9) and (10) and equations (ll) and 0-2) to possess 

similarity solutions under assumption A is that h l = “7~T 8X1,1 

h , _ where p(x L ) and q'(x 2 ) are not identically zero. 

2 q 2 (x 2 ) x 

(indices may also be interchanged uniformally in the previous statements.) 

The form of k x and k 2 indicates that the coordinate system 
(xt,x 2 ) is a modified polar coordinate system when referred to a plane. 
That is, the system ( Xl ,x 2 ) differs from the usual polar coordinate sys- 
tem only in scale variations* 

The analysis of similarity solutions for polar-type coordinate sys- 
tems is completely documented in reference 14. The coordinate system of 
this reference can be defined by h L = x 2 and h 2 = 1. The analysis 

presented in reference 14 shows that only one form for - C 2 U 2 is 

possible. This form is 

„ n mx-i 

U x = ax 2 e x 

This case is investigated further by introducing a scale change 
transformation 

X 2 

X 1 - x 2 “ e 

The equation for (ds) 2 on the surface then becomes 

2X -i 

(ds) 2 = x 2 (dx 1 ) 2 + (dXg) 2 - e 2 [(dX 1 ) 2 + (dX 2 ) 2 J 

Finally, "by an orthogonal transformation, 

X^_ = sin a Xj_ - cos a X 2 (a a constant) 

X 2 =? cos a X x + sin a X g 

The equation for (ds) 2 can be written 

C ds) 2 .e 2(Sl “ a5l + 008 a " 2 ) K ) 2 + ^ 2 ] 
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This is basically the same expression as 'gguation (36). The equation 
for U-|_ = cgUg becomes 

, 4 n'X-, + m'Xp 

= (const. )e _ 

which corresponds to equation (37). Therefore, the following important 
theorem can be stated. 

Theorem 3 . - The boundary- layer flow problems that can be solved 
by employing the coordinate system and associated mainstream flows of 
Theorem 1 are identical to those that can be solved by employing the 
coordinate system and associated mainstream flows of Theorem 2. 


Reduction of O.D.E. Conditions When Both 


Coordinate Curvatures Vanish 

The final case to be considered is k]_ = 0 and kg = 0. The co- 

ordinate lines will therefore be geodesics of the surface. In the plane 
such a system simply becomes a rectangular, Cartesian, coordinate system. 
This particular problem is discussed in reference 13 _§nd, hence, details 
will be omitted. The results of this analysis with h-j_ = hg = 1 give 

the following principal forms for U^, Ug, and (note: Ug here 

corresponds to an actual velocity component) : 


nx, i 
\ TT 1 m-1 

(l} U-^ =» ae Xg 

nxi ™ 

U 2 * be -Scf 

(2) Uq_ * ax^Xg - " 1 - 

U 2 = bxj-^ 

(3) Uq =» axq 
U 2 =» bxj 

(4) U x - ae^ 1 


= const. U]_ 


. U 1 
W const, — 

X 1 


Vl 

=s const. — 


- X 1 


» const, 
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ANALYSIS OF POSSIBLE SIMILARITY SOLUTIONS 
FOR U-,/U 2 = CONSTANT 

Throughout the previous main section it has been assumed that the 
Gaussian curvature of the surface K was identically zero. _As a con- 
sequence of this definition, it was determined that in every 

case where one of the coordinate curvatures was nonzero. In the present 
section it will he assumed at the outset that = c^Ug, and no assump- 
tions regarding the nature of K will he made. 


l 

IS! 

O 


Reduction of O.D.E. Conditions for k-j_ ^ 0, kg ^ 0 

If it is assumed that k-^ 0 and k^ 0 under the basic assump- 

tion that ~ c-jUg, the o.d.e. condition for © and d> yields 


k l = c 2 k 2 


Two possibilities now exist. One possibility is that k^ and kg 

are constant. (This was not allowed in the previous analysis except 
for kj_ ~ kg = 0.) The other possibility is that k^_ and kg are non- 
constant. If k^ and kg are assumed nonconstant,, the analysis from 

the beginning of appendix B to equation (B13) will apply since the 
analysis used only o.d.e. conditions and the results of lemmas 1 and 2. 
No recourse was made to K = 0 in that section. Hence,, by substituting 
the expressions from equations (Bl), (B4), and (B13) into equation (16), 
the following expression for K is obtained: 


K = - (kf + kl)(l + 

K - -k|(c| + 1) (l + (38) 


On the other hand, if k-i 
(16) : 


Cgkg = const., there results from equation 
K = -k 2 ( c | + 1) 


( 39 ) 



26 


MCA TM 1437 

Permissible forms for hj_, hg, k}_, kg/ ancl Ug.- - In attempting to 

determine permissible forms, we will initially distinguish between the 
case where k]_ and kg '.are- nonconstant and the case where both are 

constant. First consider the nonconstant case. 

From equations (B7) and (B13), the following relation occurs be- 
tween h^_ and hg: 



This equation can be rewritten as 

h 2 = h 1 f 10 (x 1 )f 11 (x 2 ) _ (40) 

As k^ = Cgkg, equation (BIO) is valid: 

dh dh_ 

^ = C 2^ ; (B1 ^ 

Substitution of equation (40) into (BIO) then gives 

bh l r hh l 1 

3^ " e 2 f ll( x 2>[ f lo( x l) + h l f io( x l^J 

which in turn can be written as 

S *1 y| ^ i y \ j") 

^x 2 1 " c 2 f ll( x 2) f lo( x l) ^ X]L 1 = c 2 f ll( x 2) f io( x l) ( 41 _) 

Equation (41) is a first-order linear partial differential equation 
that can be solved by classical methods. In particular, if the method 
of Lagrange is applied, a general solution for h^ is 

^l = c io f i2( x i) e<p(u) ( 42 ) 

where <p(u) Is an arbitrary function of - 

u = f 12 (x 1 ) + Cgf^Cxg) + const. 
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* 


The functions f -j g(x 1 ) and f are related to fio( x l) and f j j (x?) 
by 


f 10( x l) = f.'Tx ') J* 0 
12 ( ‘ 1 ; 

fll(x 2 ) = fi3( x 2) t 0 


The corresponding expression for hg can be determined from equations 
(40) and (42) and is given by 


h 2 = c l(fi3( x 2) e 


(P(u) 


(43) 


From equations (42) and (43) the following expressions for k^ and kg 
can be obtained: 

c. 


k-j = — <p’e -<P 
1 c 10 


k 2 = ^ < P ,e " <P 
c 10 


(44) 


If k-|_ and. kg are now considered constant, 

d In hi 


dxc 


■ h 2 k l 


Hence , 


Similarly, 


Therefore, 



In hj_ dhg 

"SZT = k l = h l h 2 k l k 2 

d^ In hg 

bx 2 m h l h 2 k l k 2 

d^ In hg d^ In h^ 
dxi dxg = dx-^ dxg 
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or - 

>.2 h 2 
o In r— 

hq 

" 0 

Trie solution to the previous equation leads directly to an equation' of 
the form of equation (40) , and the analysis that follows is applicable. 

However, with constant k^ and k 2 , the solutions of equations (44) 

can be written at once as 


<P = In (u • const. + const.) 


-1 


By means of the o.d.e. conditions it is now possible to determine 
all possible forms of the function <f> and the remaining unknowns. Be- 
fore the conditions are applied, however, the numbered coefficients in 
set (14) will be replaced by sin equivalent set obtained from the original 

uf 

list after substituting U-j_ = c£U 2 , k-j_ = e 2 k 2 , and g_ = (const.) =- k^, 
(T) and (9) , into the various terms,. This new set is 


following from 
the following: 


& 


(D* 


1 

a 

mu 2 

h l 


a Xl 

1 

a 

la (h|/U 2 k 2 ) 

hi 


dx^ 

_L_ 

a 

mu 2 hi 

h 2 


a Xg 

_1_ 

a 

in (hi/U 2 k 2 ) 

h 2 


dx 2 

1 

a 

In h 2 

*1 


K 

1 

a 

mu 2 

*2 


^x 2 

_1_ 

a 

in U 2 h 2 

hi 


axi 

1 

a 

In h-|_ 

h 2 


ax 2 


(45) 
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Proportionality between any two of these expressions will constitute 
an o.d.e. condition. 

As in the past, it 'will be necessary to make certain assumptions 
regarding the nature of U 2 before a unique set of o.d.e. conditions 

again distinguished. (Note: 


s 0 

£ 0 

^ 0 

is is presented in appendix E. The re- 
sults can be summarized in the following theorem. 

Theorem 4 . - Let (x^,x 2 ) be an orthogonal curvilinear coordinate 
system for which kj_ ^ 0 and kg 0. Further, assume that 
U-j_ = (const. )U 2 . Then, necessary and sufficient conditions for equa- 
tions (9) and (10) and equations (ll) and (12) to be reduced to ordi- 
nary differential equations under assumption A and the given similarity 
transformation are that 



where u = p(x]_) + c 2 q(x 2 ) + c 3 . (The symbols c-j_, c 2 , and c 3 
not to be confused with previous cases where these constants were 
employed.) 

(2) The following sets of conditions hold for U 2 and <P : 
(a) U 2 = constant: <P = ln(au + b) n or <p = au + b 


can be obtained. Four cases aa 
Up = c L U 2 ) . 

(1) Case A: U 2 = const. 

Su 2 , SUo 

(2) Case B: _ fi 0; 






(3) Case C: s 0; 


2 u ~ 2 


(4) Case D: ^ ? 0) 


The analysis of these four ease 


are 
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,, = s< l( x 2) 

v c ) U 2 = re : <P = au + b 

(d) U 2 = re s<P : <P = au + b or <p = ln(au + b) n 

_ sp( Xl )+t<i(x 2 ) 

(e) U 2 = re : <P - au + b 

( 3 ) S = (const. 

where s and. t are constants • 


,, R+1 ^ + t J ±S , p ?i nt ' it is of Merest to determine the cases that are 
it lvi ?° 3e refer £ ed to ln Theorem 1. At the outset, note that 

’ diere * ■ au + b *» «* £s - s* 


a Pl x lJ_ac 2 q(x 2 ; 
1 e 


hi = (const. )p'(x-L)e 

- (const, ) e ac 2l( x 2) A[>( x l>] 

h 2 - (const. )q.' (x 2 )e P ^ x l^ e ac 2 q t x 2) 

- (const.)e ap(x l ) ^[ e ac 2 l( x 2 >] 

By denoting e aP(Xl) by P^) and by , ( ), ^ 

expressions can be written as 

h l = (const. )P , ( x 1 )Q(x 2 ) 

h 2 = (const. )Q’( x 2 )P(x 1 ) 

OJhese expressions are the same as those used to define h, and h in 

^ltiS iram X \ f ° 11 2 W8 ^ th ® re f ore ^ that K = 0, and the solutions re- 
sulting from <p - au + b duplicate previous solutions. However it is' 
readily verified that solutions different from those given in Theorem 1 

result when <P = In (au + b) n . For n/o,K / 0, the correspond!™ 
th^orm! arS n ° ndeVel0pable - 111686 observations lead to. the fallowing 
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Theorem 5 . - The forms for h-^, hg, and Ug given in Theorem 4 
which are associated with nondevelopable surfaces are characterized by 

cp = ln(au + h) n . The expressions for h l* h 2> and Ug can "be written 

h^_ = (const. )p r (x-j_) [• •p( Xl ) + b' 4 ( 1 , 5 ) + =’] n 

hg = (const. )q.' (xg) + b'q(xg) + c'J n n^ 0 

U 2 = (const. ) ^a'p(xjL) + b'q(xg)J m 

Geometric considerations. - The specific expressions for h l> h 2 ' 
and Ug given in Theorem 5 make possible a rather general analysis of 
the geometric aspects of the problem. In order to simplify the analysis, 
P( x l)j the constants that appear in the expression for h-^ 

and h 2 are chosen in such a vay that 

h-L = (aX]_ + ^ x 2) n 

h 2 = (bx-l + bxg) 11 

Given the corresponding differential quadratic form, 

(4 ») 2 = (a^ + b *,) 211 ^) 2 + (dx 2 ) 2 ] (46) 

the existence of an associated class of surfaces is assured (ref. 21 , 
p. 122) . This class of surfaces will now "be studied more closely. Con- 
sider an orthogonal transformation of coordinates defined by 


= X-^ sin a - Xg cos a 
¥: 

= X^ cos a + Xg sin a 


(a a constant) 


Under this transformation, equation (46) becomes 

(ds ) 2 = £(a sin a, + b cos a)x£ + (b sin a - a cos a,)X^j n £(dX *) 2 + (dX^) 2 ! 


If a is now chosen such that 


a sin a + b cos a = 0 


there results 


(ds ) 2 = c 2 (X*) 2n [(dX *) 2 + (dX*) 2 J 
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where 


Finally, "by choosing 


c = h sin a - a cos a 


x 2 = 5 Hhr < x 2) n+1 ( n ^ - 1 ) 


n 


xx - o(^r x* 


equation (47) can be written 

l 2 _ v 2m/,„ \2 


(ds)^ * Xg (dX^^ + (dX 2 )' 


where 


m = 


(48) 


n + 1 

It follows frcm the above equation that m_ . cannot be equal to 1. 
If n = -1, 

X 2 = c' In X| 

Xj.= c'Xj 

The differential quadratic form (47) becomes 

(ds) 2 = e " 2X 2/ c ’ (dX x ) 2 + (dX 2 ) 2 


(49) 


Now consider surfaces in space in which it is possible to embed 
coordinate systems having differential forms corresponding to equations 
(48) and (49). First note that hq in equations (48) and (49) is a 
function of X 2 alone. Now for a (ds) 2 of the type 

(ds) 2 = f 2 (x 2 )(dXq ) 2 + (ax 2 ) 2 

a specific class of surfaces can be associated (see ref. 21, p. 206). 
These surfaces are surfaces of revolution given by (Yq Cartesian 
coordinates) 


q; r?> 
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X-, 

Y-l = af (Xr,) cos ~ 

Y 2 = af (X 2 )sin -i 
8 . 

Y 3 = f a/i - (a) 2 (f') 2 dX 2 + const. 


Ia relation to the present problem this report w ill first investigate 
surfaces for which f(X^) = Choosing = 0 and eliminating the 

remaining parameter show that the surface of revolution is generated by 
rotating the curve: 


- 1 dY. 


(50) 



about the Y^-axis. This is also the class of surfaces investigated "by 
Geis in reference 7. 


The nature of the curves can "be qualitatively determined from an 
examination of the derivative 



For 0 < m < 1 the curve should have the following general shape: 




34 




For this particular case, it is known (ref. 21, p. 207) that the surface 
of revolution is obtained from revolving the tractrix shown below about 
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n 

i • 

.•a 


As indicated in reference 7, the characteristics of the surface are such 
that the "boundary-layer equations will prohahly apply over lim ited regions. 
It is, of course, necessary to keep in mind that, in any flow problem 
over curved surfaces, the "boundary- layer equations are valid only in re- 
gions where the minimum value of the principal radii of curvature does 
not exceed the "boundary- layer thickness. 

While a special class of surfaces has "been exhibited with coordinate 
systems satisfying the expressions for (ds ) 2 in equations (48) and (49), 
it is well to consider what other surfaces might be admissible. The 
admissible surfaces are surfaces that are "applicable” to the surfaces 
of revolution illustrated previously (ref. 21, pp. 172-174). This means 
the class of surfaces obtained from. the surfaces of revolution by bending, 
without stretching, compression, or tearing. 

The streamlines of the mainstream flows will be lines that make con- 
stant angles with the coordinate lines. In systems where one set of 
coordinate lines corresponds to meridian curves on a surface of revolu- 
tion, the stream-lines are the loxodr ernes of the surface. 


Reduction of O.D.E. Conditions When One Coordinate 
Curvature Is Identically Zero 

Solutions of the o. d. e. conditions are now to be obtained under the 
assumptions that one of the coordinate curvatures vanishes identically 
and that = c-^Ug. 

Initially, consider the case of k^_ = 0. (The case for k 2 = 0 then 

can be readily determined from symmetry considerations . ) Once again, a 
set of o.d. e. conditions is obtained by assuming various possible forms 
for U 2 . 

The analysis is presented in appendix F and can be summarized in the 
following theorem, which is stated for k^_ 0 and k 2 = 0 . 

Theorem 6 . - Let (x^,Xg) be an orthogonal coordinate system for 
which k-j_ 

necessary and sufficient conditions for equations ( 9 ) and ( 10 ) and (ll) 
and ( 12 ) to be reduced to ordinary differential equations under assump- 
tion A are 

( 1 ) U 2 = const., and 


^ 0 and kg = 0 . Further assume that TJ-i =■ (const. )U ? . Then 
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or 


(2) 


h l = 

<3.' £x ? ) 

h2 -«C& 


h l = 

_ <l( x z) . , 

h 2 " n ^ere a f 0 


U 2 . (=o n st.)e^ P(Xl,aXl 

h - P(X 1 } 

1 " <l( x 2 ) 

q' U 2 ) 


(3) Ug a* (const. Jq^Xg), and. h-j_ and hg are expressed "by either 
the first or second set of relations in condition (l) 

(4) Ug = (const. ) q n (xg,)e */'p( x i) ( ^ x 1^ and and hg are expressed 

as the set of relations in condition (2) 

The function g2 is determined by 

g 2 = ‘(const. )Ugk^ 

( Indices may be unif ormally interchanged in the previous statements. ) 

Geometric considerations . - If the Gaussian curvature is calculated 
from the values of h^ and hg in condition (2) of Theorem 6, it can be 
shown that K = 0. Consequently, cases (2) and (4) are completely 
covered by previous analyses . 
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in 

e'- 

en 


The first set of relations for h-^ and h 2 in condition (l) of 
Theorem 6 can he written 


h^ = e 


nxg 


h 2 = ! 

hy choosing p^x-^ = 1, q-^Xg) = e *2. 

The quadratic differential form then hecames 
(ds) 2 = e^SCd^) 2 + (dx 2 ) 2 


This form is identical to equation (49) with c’ = - jj. 

The second set of relations for h^_ and h 2 in condition (l) of 
Theorem 6 can he written 

h i= 
h 2 = i 

hy choosing 


qi(x 2 ) = -nx" 1 
Pl ( Xl ) = (-n)- n 


For this case. 


( de) 2 = ( Xj)- 211 ^) 2 + ( to 2 ) 2 


This result can he compared with equation (48) with m = -n. However, in 
this equation it is possible to choose n = -1, which leads to a devel- 
opable surface. (The value m = 1 was not allowed in eq. (49).) Refer- 
ring to equation (50) shows that, for m = 1 and — g- = 1, ~ ~ 

and the corresponding surface of revolution is a cone > lj (i 
lowing sketch) or a circular cylinder = 1^ . 


er 

'see fol- 
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As would "be expected, the expressions for Uj_ = c^Ug determined from 
Theorem 6 can readily be shown to be identical to expressions for 
determined frcrn Theorem 5 after the form for (ds) ^ in Theorem 5 is 
transformed into the form determined by Theorem 6. 


This section concludes by observing that further analysis consider- 
ing = 0 and k 2 = 0 is not necessary. In this instance K = 0, 
and this case has been completely analyzed. 


GENERAL SUMMARY OF REQUIREMENTS FOR THE EXISTENCE OF 
SIMILARITY SOLUTIONS 

The various results obtained for flows over developable surfaces 
and for flovs_over nondevelopable surfaces characterized by 
Ui = (const. )U 2 can be summarized as follows. For surfaces of noncon- 
stant Gaussian curvature, constant curvature, and zero curvature, a 

suitable choice of Gaussian coordinates can be made in each case that 

reduces the corresponding forms of h-j_, h 2 , Uj_, Ug, and to a basic 

form given later. Associated with each basic form there can be .chosen a 
characteristic surface of revolution. The meridian curves on the surface 
and their orthogonal trajectories constitute the coordinate lines of the 
basic system. All other permissible coordinate systems are obtained from 
the basic system by a scale transformation and an orthogonal transforma- 
tion. All permissible flow surfaces for a given value of K are the 
class of surfaces applicable to the characteristic surface of revolution. 

The basic forms for h^, hg, U^, Ug, gZ, and the equation for the 
associated surface of revolution are as follows: 


(1) K / constant 

hi = x^ n ^ 0,1 
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h 2 = i 

U-l = (const. )Ug = Xg 

2 U-, 

g = ( const. ) — 

*2 

Surface of revolution: obtained by revolving the curve 



- 1 cLY^ + const. 


about the Y 3 -axis. 

(2) K = constant (nonzero) 


h L = e^ 

h 2 - 1 


mxg 


U-j_ = (const. )Ug = e 
g 2 = (const. )U]_ 

Surface of revolution: obtained by revolving the tractrix 


Y 3 = a - 1 ail + const. 


about the Y 3 ~axis. 
(3) K = 0 

(a) kf + k| 4 0 


* 2-1 


_ mx-i 

U x » (const. )Ug = Xge 

g 2 * (const. ) 

x 2 


Surface of revolution: a cone 
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(*) ^ 


k 2 = 0 


h 1 = h 2 * 1 


1. U, 


nx-, m-1 
5 lx 2 


TT , nxi 
U 2 = be 

g 2 - (const. 


2. U, = 


n m-1 


tt -u a-1 m 
U g = bx-L xg 


g 2 = (const. ) 


Ul 

xi 


3. U n = 


n 

ax-^ 


U 2 = bx* 


U, 


g 2 = (const. ) -i 


4. U-^ = ae 


nx]_ 


U 2 = be 


mx^ 


g 2 as (const. )U^ 

Surface of revolution: circular cylinder 


METHOD AMD RESULTS OP REFERENCE 16 

The assumption of proportionality among coefficients in equations (9) 
and (10) and (ll) and (12) (assumption A) was used in preceding sections 
as a sufficient condition for reducing these equations to ordinary dif- 
ferential equations. Geis (ref. 16) has shown in the meantime that such 
proportionality is also a necessary condition for the complete boundary 
value problem. He actually shows that the previously mentioned coeffi- 
cients divided by g 2 are constant. The proof is based on a detailed 
study of these coefficients when the basic equations are subject to the 
boundary conditions on F( T| ) and G(t|). This proof is the content of 
section 4 of Geis' paper. 


SD6? 
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The search after permissible forms of the functions h.^, h 2 and 
Ux, U 2 proceeds in the present report along lines that are quite differ- 
ent from what is done in the corresponding section 5 of reference 16. 

Geis discusses a subcase of the case K = 0, namely, = kg = 0. The 
rest of his discussion uses the alternate assumption that 


d(h,U,) -s 
' Sxg' - (h 2 U 2 ) = 


which is satisfied if the main flow is irrotational. This hypothesis 
stands in contradistinction to the assumptions of developable surfaces 
or = (const. )U 2 made here. It is remarkable that all but two of 
the solutions obtained by Geis for U^, Ug, h^, h 2 , and g 2 are also 

found in the present analysis. In two exceptional cases Geis gives 
implicit solutions that were not obtained here. The first of these 
solutions is characterized by 


*1- 1 


h X - l 


TT 

2 " V 1 - V 

h 2 = =- 
2 Uo 


g 2 = i- 
& 2cp 


where cp = cp(x-j_) is a solution different from (x-|_ + const.) of 

<pcp" + (<P' “ l) 2 = 0 

The second of these solutions is characterized by 

U 1 = <P( x l) U 2 = A/^uJ + B 

g 2 - -cp' 

hi=l h 2 = ^ 


where cp is a nonconstant solution of 

(cp 2 + B)cp" + 0T <l> 3 " Cep 2 + \ Bcp - BC^cp’ 2 = 0 

Geis notices, however, that his explicit solutions are also valid when 
his condition of irrotationality is not satisfied. 

It might finally be added that certain other solutions given im- 
plicitly in reference 16 are found in explicit form in the present paper. 



42 


MCA TM 1437 


It mi gh t he mentioned again that the .emphasis here has been to 
examine the geometric aspects of the problem, that is, to determine the 
nature of coordinate systems, main-flow streamlines, and flow surfaces 
from the solutions found for Uj_, U 2 , h^, and hg. This would appear to 
be a necessary requirement for fully evaluating the physical significance 
of the results. 

UNCOUPLING OF EQUATIONS 

The principal advantage of employing similarity techniques in an 
analysis of boundary- layer flows is the reduction of. the partial differ- 
ential equations for the flow to ordinary differential equations. How- 
ever, the solution of the system of ordinary differential equations is 
generally difficult to obtain. Inspection of equations (9) and (10) or 
equations (ll) and (12) discloses that the equations "are of order three 
and nonlinear. Solution of systems in the past has visually required 
application of numerical techniques. If a wide range of flows is to be 
studied, this can be laborious and time-consuming. An additional diffi- 
culty that arises is that in a large number of cases the equations are 
coupled, which means that the functions F(q) and g'(tj) appear in both 
equations. If the equations are uncoupled, and one equation is expressi- 
ble in terms of a single function, that equation canTse solved quite 
readily, and the values obtained can be used in the solution of the re- 
maining equation. It is of interest to know, therefore, for what classes 
of flow problems this is possible. 


With this concept in mind the present section will determine the re- 
strictions on hj_, hg, U-p and U 2 that lead to the uncoupling of equa- 
tions (9) and (10) and (ll) and (12). 

Uncoupled equations in the trivial sense of one equation identically 
vanishing will not be considered. (This could cane about in certain . 
classes of flow problems in which = 0 or Ug = 0 and the main-flow 
streamlines have zero geodesic curvature. ) 


The requirement for uncoupling is simply that the coefficients of 
terms involving both F and G should vanish in either equation (9) or 

(10) (or eqs. (ll) and (12)). Without loss of generality the requirement 
for uncoupling will be imposed on equation (9) (results hold also for 

(11) ). Examination of equations. (9) and (ll) discloses that the follow- 
ing equations must be satisfied: 


Up 5 In U-i 

H ~ + ^ S ° 


( 51 ) 
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-i 3TJ, Up ^ 2 

— r— - - ■ . ^ n - 5 ■ + Upkn 

hg dxg 2hg oxg 2 1 


0 


*1 


_ kp = 0 


Equation (53) gives at once the following lemma. 


(52) 

(53) 


L emma 5 

The ordinary differential equations (9) and (10) (or (ll) and (12)) 
are uncoupled only if one of the coordinate curvatures vanishes identi- 
cally, that is, only if one set of coordinate lines is a system, of 
geodesics on the flow surface. 

Attention is first restricted to flows over developable surfaces. 
First, note from equation (51) that, if a In Uq/dxg = 0, then k-j_ = 0 
and, hence, both sets of coordinate lines are geodesics. If it is now 
assumed that a In U^/Sxg ^ 0 and k]_ 0, then ( see appendix D) 


Uq = Cl u 2 

Hence, from o. d. e. conditions and CD in set ( 14) , 

g 2 = CgU-,]^ 

. a In g 2 = * ^ U 1 + 5 ^ k i 

*' Sxg Sxg dxg 

But from equations (51), (52), and = c^Ug, 

d S 2 = o 
Sx 2 

a In Ui a In ki 
axg ~ SET 

Examination of equation (17) shows that for K = 0 

2 a In k-L 

*2 Sx 2 = " k l 

Substituting this result in equation (54) gives 

1 5 111 °L . 

hj 5^ ‘ 


(54) 
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If this expression is in turn substituted in equation (51) , 

2*1 = 0 __ . 

which contradicts the assumption that ki ?= 0. The following lemma now- 
results. 


Lemma 6 

For flow over developable surfaces, equations (9) and (10) are un- 
coupled only if ]£}_ s kg = that is, both sets of coordinate lines are 
geodesics. 

If k]_s kg = 0, it follows at once from equation (51) that 
U 1 = U^(x-j_). Furthermore, examination of possible forms for U^, Ug, and.^. 
g2 for this case (see p. 40) shows that in every instance g2 = (const. ) t^* 
or g2 = (const. Hence, with Ui = ^(x^) equation (52) becomes (now 

it is required that Ug = Ug) 


Su 2 


= 0 


-J 

cn 


Therefore, the permissible flows can he reduced to the following two 
cases (see again p. 40) : 


nx-i 

U x = ae 1 
U 2 =be 1 


or 


U-^ = ax^ 


U 0 


bx 


m 


Now, consider flows over nondevelopable surfaces with K const. 
As kg s 0, the permissible solutions for hi, hg, g 2 , and Ui can be 

taken as those given in case (l) (p. 38) : 

*1 = *2 1 

ho =1 r_ ■ - ■ 


^ = (const. )U 2 = x“ 
g 2 = (const. JU-^/xg 
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Now, for Uq_ = (const. )Ug it was previously shown that 

5 ^ S 2 . = 0 

dXg 

However, = (const. and, therefore, it is necessary that m = 1. 

Substituting the permissible forms for U-. , h-, , and h ? into equation 
(51) gives 


JL_ 

*2 



0 


Hence, n = -1. 


While the previous . results appear to lead to a valid case of un- 
coupled equations, substitution of the actual forms into equation (9) 
shows that this is not so. As may be verified, substitution leads to 
the equation 


F ,n (tj) = 0 

The function F( t]) cannot, therefore, fulfill the imposed boundary con- 
ditions, and this case must be ruled out. 

Finally, consider K = constant (nonzero) . The solutions for *1> 
hg, g^, and U^_ of case (2) (p. 39) then apply. 

As 

g2 = (const. ) U-j- 

and, as once again 

5 jif S 2 = 0 

SXg 

there results 

S In U-, 

= 0 

This result substituted in equation (51) yields kq = 0. This, however, 
would require (with kg = 0) that the surface be developable and would 
therefore contradict the hypothesis. 

The results are summarized in the following theorem. 
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Theorem 7 

The ordinary differential equations (9) and (10) are uncoupled if 
and only if 

(l) Flow surface is developable, and both sets of coordinate lines 


are geodesics. 




(2) U-]_, U 2 , and g 2 have one of the ’two following forms (except for 
change in indices) : 

nx - - 

I U = ae x 

& 

Si 

(A) 

1 mx-L . 

U 2 = be 


■ — 


i g 2 = eU 1 “ 

- 

• — 


' 1 J 

1 

U 

11 


. 

(») ! 

°2 - K 

. 



2 U 1 “ 
1 8 - c 7 

X 1 

- ■■ ■ 


Equations (9) and (10) have, respectively, the following forms if (A) is 
employed: 



nj(F') 2 _ . ij - cF" 1 = 0 

(55) 



m(F'G' - 1) - nG " F - cG 1 " = 0 
2 

(56) 


If the forms of (B) axe employed, there result, respectively, 
tions (9) and (10) 

for equa- 



n |(F' ) 2 - lj - (n + l) - cF" 1 = 0 

(57) 



m(F’G’ - 1) - (n + 1) - cG"’ = 0 

2 

(58) 


For both sets of equations the boundary conditions are 


w 


F(0) = F'(0) = G(0) = G' (oy = 0 


stm 
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11m F’(tj) = 1 11m Q'(ri) = 1 

Tj -► oo T] 00 

It is not necessary to consider equations (ll) and (12) as Ug = 0 
leads to straight-line flows in the present case. 


PRACTICAL APPLICATIONS OP THEORY 

Pram, a practical standpoint, similarity solutions might he applied 
to the study of boundary- layer flows over such aerodynamic configurations 
as wings, missiles, fuselage forms, or channel flows. Pram the stand- 
point of analyzing flows over wings or in channels, the only type of 
analysis that seems promising is the one employing rectangular coordinates 
(i. e. , coordinate lines are geodesics). The principal reason is that in 
such configurations the boundary layer is generally initiated along a 
straight line on the surface (e. g. , leading edge). This physical case 
can only be approximated when an analysis allows a boundary layer of zero 
thickness to exist along such a line. This obviously can only come about 
in a rectangular coordinate system analysis. On the other hand, polar 
and spiral coordinate systems’"may have particular application to flow over 
such configurations as missiles where the boundary layer develops from a 
point (e. g. , nose of the missile). 

In any of these cases, however. It should be kept in mind that a 
similarity analysis will generally predict only qualitative behavior of 
the flow. The restrictions imposed on the main-flow velocity components 
will, in general, be too severe to conform to a specified flow configura- 
tion, and, at best, an approximation to this flow can be constructed. 
Nevertheless, experimental verification of certain aspects of flow behav- 
ior predicted by theory has been very encouraging in at least one in- 
stance. The investigation presented in reference 1 shows that limiting- 
flow deflection on a channel surface can be accurately predicted. 

The calculation of boundary- layer velocity profiles from a similarity 
analysis can also serve as a guide in setting up approximate analysis of 
boundary- layer flows using so -call led mcmentum-integral methods. This 
method generally requires an a priori specification of the velocity pro- 
file shapes, which are then approximated by an analytic expression. Lack 
of information on the forms of three-dimensional velocity profiles has 
seriously hampered application of this technique. 


CONCLUDING REMARKS 


The following conclusions can be drawn from the analysis presented: 
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1. The requirements for a similarity analysis of the boundary- layer 
equations for flow over developable surfaces can he completely determined 
if the equations are referred to orthogonal coordinates. Two basic sets 
of solutions are sufficient for analyzing all permissible flows. 

2. For mainstream flotfs over nondevelopable surfaces with U^_ = cjJJg, 
two basic sets of solutions are sufficient for all permissible flows. One 
set applies to surfaces with K ^ const., while the other applies to 

K = nonzero const. All permissible flow surfaces are surfaces geometri- 
cally applicable to two basic classes of surfaces of revolution. 

3. Uncoupled systems of ordinary differential equations resulting 
from a similarity analysis occur only when coordinate lines are geodesics 
( rectangular system in the plane). 

In conclusion, the following problems pertaining to the analysis of 
similarity solutions seem worthwhile for continued investigation: 

1. A general analysis for all possible types of similarity solutions 

has not yet been evolved in the sense that investigations to date (reft 
16 and the present paper) have employed certain assumptions to reduce the 
complexity of the problem. _ __ 

2. Solutions of the ordinary differential equations arising in a 
similarity analysis are few in number. Little work has been done on 
coupled equations. It would be of interest to study variations in solu- 
tions for a range of parameter values either through a program on high- 
speed computing equipment or by developing .suitable approximation 
techniques. 

3. Extensions of the present theory for incompressible laminar- 
boundary-layer flow to compressible laminar-boundary-layer flow would 
be of interest. Seme investigations of this kind have been carried out 
(e.g. , ref. 22). 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, June 10, 1958 
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Ai 

a y a ^ 

\),Td ± 

c i 

d,d ± 

e 

F, P(ti) 

f ' f i 

G, G(ti) 

Sij 

gjgCx^Xg) 


-u* 

h l’ h 2 

K 


k l> k 2 


m 

n 


APPENDIX A 
SYMBOLS 

constants 

constants 

constants 

constants 

constants 

constants 

constants 

constant 

functions of tj 

arbitrary functions 

functions of T) 

metric tensor for three-dimensional coordinate system 

arbitrary function occurring as a factor of tj 

square roots of metric-tensor components in orthogonal 
coordinate system 

specific forms for hj_ and hg defined by eq. (29) 
Gaussian curvature of surface 

geodesic curvature of coordinate lines of x^_ and xgj 
respectively 

constant 

constant 
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Pi 

<J± 


s 

t 


Ui,u 2 

u? 


U 2 


u 


“i 

V 


Xi,X 2 
v * Y * 

Xi,x 2 

Y 1' Y 2' Y 3 

y 

y* 

n 

e,p 

V 

Si 

cp(u) 


arbitrary functions — 

arbitrary functions 

constant ~ 

arc length 

constant 

mainstream velocity components in x-j_- and x 2 -directions 

function of x^xg; eq. (6) 

common designation for either U 2 or U 2 

special designation for function of x^ and x 2 occurring 
in eq. (42) 

boundary- layer velocity components in x^- direct! on 

transformed boundary- layer velocity component normal to 
surface — 

coordinates defined by eq. (27) 

curvilinear coordinates ~ 

coordinates 

curvilinear coordinates 
Cartesian coordinates 
y*/Vv 

physical coordinate normal to flow surface 

similarity parameter _ ” 

coordinates defined by eq. (C5) 

coefficient of kinematic viscosity 

function of X]_,X 2 , eq. (Cl) 

arbitrary function, eq. (42) 
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Subscripts : 

i, j,k,Z,m,,r,6 denote index numbers 

Superscripts : 

Primes denote differentiation 
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APPENDIX B 


ANALYSIS : K = 0; fk 0, kg jk 0. PROOF OF THEOREM 1 

The solutions for the various unknowns in set (l4) on page 10 are 
obtained here for the case K = 0; k^_ 0, kg ^ 0. 

Permissible Forms for k^_ , kg, h^, and hg 

Before k l' kg, h-L, and hg can be determined, certain other 

general results are needed. In this regard, it wild first be shown that, 
under the assumption of neither kj_ nor kg being identically zero, it 

must follow that k-)_ and kg each possess nonvanishing first derivatives 
with respect to both x-j_ and xg. 

Assume that a derivative of either k^ or kg 'with respect to one 
coordinate is identically equal to zero, while the derivative with re- 
spect to the other coordinate is not. 

Ski Ski 

Specifically, assume t- — g 0 and — — s 0. It will first be shown 
' 8k, Sx l 

that the assumption ^ — ' U 0 leads to the equation ■ ■ 

Sx 2 . , d. ' 

2 (B1) 

h 2 Sxg 

From lemmas 1 and 2 and equation (19), 

2 c|. 

g = — Ugk-L 


Hence, 


Sing 2 ^ 1x1 U 2 ^. lnk l 

= < r - + — <r— 


Now, both S In g 2 /Sxg and S In TJg/Sxg cannot be zero, for then 

S In k-7 

— ^ = 0, which would violate the initial assumption. However, if 

2 

S In S In Ug 

— g - - ** ■ = 0 and — jk 0, from equation (B2), o.-d.e. condition for 

Xg Xg 


316 ? 
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(6) and (?) , and lemma 1 there results 

d In k^_ d In U 2 

dx 2 ~ ” dx 2 - ^l^^l 

It therefore follows that equation (Bl) holds. 

5 In a 2 S In U 2 

If — ^ 0 and — = 0, then from equation (B2), o.d.e. 

condition for © and <D , and lemmas 1 and 2, 


5 111 k l 5 In g2 

■v = ■ - ■ ‘ V ' ' = 


dx 2 Sx 2 

and, again, equation (Bl) holds. 


t 2 il 2 k l 


Finally, if both d In ^ fbyi^ and ^ In Ug/Sxg are not identically 

zero, equation (B2), the o.d.e. conditions for (5), (6), and (7), and 
lemmas 1 and 2 yield 


5 In 


ox? 


- 


and equation (Bl) holds. 

Now, from equation (Bl) and the assumption that 

S 2 In ki 


cik-, 

m °' 


Bhg u jji ±s.-j^ 

k i^= Coonst ^ ax ~ ax 2 ~ E 0 


She 


However, k-j_ ^ — = h-^tLgk^kg = 0 contradicts the basic assumption that 

neither k-^ nor kg vanish identically. Similar ly, it can be shown 

Ski a Ski 

that the case ^ — p 0 and ^ — e 0 is incompatible with that 
1 "^2 

assumption. 

S imi lar statements hold for kg. Hence , the following lemma can 
now be~ stated. 

Lemma 5 . - If both k^_ and kg are different from zero and neither 
reduces to a constant , then both coordinate curvatures kj_ and kg must 
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possess nonvanishing first partial derivatives with. respect to both x^_ 
and x 2 . 

By employing a procedure similar to 'the one used in establishing 
equation (Bl) , it is also possible to show that 


or, by using lemma 2, 


dg 3 In k-j_ 

kl = ^1 dx l 


d£ d In kg 
kz = ^ X 1 


Hence, from equation (B4) and the definition of kg, 

dg 3 In kg 3 In hg 
c 2 ~ ^1^2) 


Therefore, 


3 

3x7 


^ ^ h 2° 2//d2 ) = 0 


••• ^2 = h 2 2/d2f l( x 2) 


where f^(xg) is arbitrary. 

From the definition of k^, equation’ (Bl) , and lemma 2, 

3 In h x 3 In k-j_ 3 IS k 2 


h 2 k l = d 


1 3xo d l 537 


(B3) 

(B4) 


(B5) 


••• ^ ^ Vi 17 - 1 - 

1 / d l / \ 

k 2 = h-L f 2 (x 1 ) 


(B6) 


where fg(x^) is arbitrary. From equations (B5) and (B6) is obtained 

rf,(x,)i 42/o2 

(B7) 


= h d 2/ d ! c 2 


f 2 (*i)T 

■l( x 2l. 


Now, relations between the various constants appearing in equation (B7) 
must be found. 


n t 
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First of all, by differentiating with respect to x-j_ 

according to equation (Bl) , 


1 L Sh 2 v. Sk i\ 1x1 k l 

%; K si + h 2 sr 


<^ x l dx 2 


(B8) 


Also, from equation (B3), 


h l k l 


5 In k]_ 
dxi 


and differentiation of this expression with respect to x % gives 

2_ ( dh^ dk-| \ d 2 In k-^ 

^ \ k i + h i si/ - S2 ( B9 ) 

From lemma 2 and the definition of k-i and kg, finally, 

din dh,. 


dx 2 C 2 dx-^ 


(BIO) 


In equation (B9) substitute for Sh^/Sxg and dk-]_/dx 2 according 
to equations (BIO) and (Bl) : 


dh P kfh? \ S 2 hi k 


i klC2 ^2 + fp h 

d 2 l 1 ^ ox^ dj_ J-y dx-^ 0X2 
Similarly, equation (B8) becomes, by equation (B3), 


(Bll) 




Sh 2 2 h x h 2 \ d 2 In ^ 


1 ^7 + k l ” 


^- 2 


(B12) 


Finally, equations (Bll) and (B12) yield 

k l c 2 ^ h 2 _ ^1 ^ h 2 
d 2 dxj_ “ dj_ Ski 


Hence, 


c 2 

djAl- 1 


(B13) 
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It is now possible 
in equation (17) may be 
(B4): 


to evaluate cU by noting that the derivatives 
eliminated by the use of equations (Bl) and 


kf + 


k 2 

*2 


+k t? 



which, by equation (B3), becomes 

(kj + + ° 

whence d^ = -1 as k-j_ and kg are assumed nonvanfshing . 
It therefore follows from equation (B13) that 

g.-i 

*2 


From the previous results, equations (B5) and (B6), and lemma 2 there 
results lemma 4. _ 

Lpmrnfl. 4 . - Under the assumption of nonyanishing coordinate curvature 
tures, and kg must be expressible as l- - 

kg = fg(x 1 )/h 1 


k l = c 2 f l 

Expressions for and hg can now be obtained as follows . From 

the definition of k-j_ and kg, 

5 In ho 
d In ho 

&c" '~ = c 2 f l( x l) ( B15 ) 

Hence, h^_ and hg must be expressible as a product of a function of 
Xj_ and a function of x 2 : 

h x = q 1 (x 1 )q 2 (x 2 ) (B16) 

hg = PqCx^JpgCxg) (B17) 
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1 


XI 

i 




Substituting expressions (B16) and (B17) into equation (BIO) gives 

= c 2 p£(x 1 )p 2 (x 2 ) 


Neither q^(xg) nor p-^(x^) can be identically zero under the assumption 
that neither of the coordinate curvatures vanishes identically. Therefore, 


<3.^X1) = c 3 p|(x 1 ) 
P 2 (x 2 ) = c^Cxg) 


Hence, h^ and hg can also be written in the form 

h l * c 3Pi( x i)<l2(x2) (BIB) 

^2 = c 4Pl( x l)^(x 2 ) (B19 ) 


From the definition of and k 2 

■i 5hp 

kp = — ±— —A = 

h^ a*[ 



and equations (B18) and (B19), 


1 


(B20) 


1 

c 4 PlCxi)q’(x2) 


(B21) 


Permissible Forms For U, , 1J 2 , and g 2 

It is noted once more that before a unique set of o.d.e. conditions 
can be prescribed for the determination of U-j_ and U 2 , certain assump- 
tions regarding the vanishing or nonvanishing of various terms in set 
(14) must be made. As the terms of interest in set (14) involve the 
first derivatives of and U 2 , the necessary assumptions can be im- 

posed in terms of the vanishing or nonvanishing of these derivatives . 

By re calling from lemma 1 that = CjUjj, the following four cases 
cover all possibilities: 

dU-, dU 2 

(1) Case A: ^ = C 1 ^ 4 0 

dU-L dU 2 

dx 2 _ C 1 Sx 2 



and 
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(2) 

Case 

B: 

kx-j_ 

SO 

bU ? -- - 

— s 0 — - - - — 

ox^ _ 




dU x 


5u 2 , 


and 


^2 

= C 1 


(3) 



dU-, 


dU 2 , 

Case 

C: 

SEJ 

= C 1 

'd^£° 




dU-L 


au 2 


and 


^2 

s 0 

55 p 0 _ ■ - 

(4) 

Case 

D: 

U 1 

and 

U 2 are nonzero constants (all derivatives 


vanish) . 

First consider case A: 

From the o.d.e. condition for © and © and lemma 1, 

Bu 

T ± 5^ - C 5 U 1*2 

a 

•'* ~ = c 5*2 h l 

Employing the expressions for h-^ and hg given in equations (B16) and 
(B20), respectively, gives - ■ - 

a In u x d In p 1 (x 1 ) 

~ c 5 dx^ 

In U x = In [Pi(x;l)] C5 + ^(* 2 ) ^ B22 ^ 

Now, from the o.d.e, condition for (£) and © and lemma 1 , 

a In 0 L 

- = e zV8 

As above , In Uj_ can be evaluated to obtain 

In U x = ln(q 2 (x 2 )] 6 + f^x^ (B23) 
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It then follows from equations (B22) and (B23) that U-^ and U 2 must 
be expressible as 


U 1 = c 7Pl( x l)l2C x 2) = c l u 2 


(B24) 


where Cj, n, and m are nonzero constants. 

The form of the function g 2 is obtained from the o.d.e. condition 
for and (6) and lemma 1: 

g 2 = Cg^kg 


Hence , 


v 


g 2 = c g pj _1 (x 1 )q“- 1 (x 2 ) (B25) 

If the expressions for k-^, kg, h^, hg, U-^, Ug, and g 2 are substituted 

into remaining o.d.e. conditions, it will be seen that these conditions 
are satisfied. Hence, for this case, necessary and sufficient conditions 
for obtaining ordinary differential equations from the original partial 
differential equations under assumption A have been found. 


The analysis for cases B, C, and D follows the method outlined 
previously. The results are as follows: 

Case B: = (constant) q 2 (x 2 ) = c-jUg 


g 2 = (constant) 


m-1, x 
q 2 (x 2 ) 


Pl(x;]_) 

Case C: U-^ = (constant) p^(x-]_) = c-]U 2 

Pl" 1 ( x l) 


g 2 = (constant) 


<l2(x 2 ) 


Case D: U-j_ = (constant) = C]_U 2 


S 2 = 


Constant 

Pl(Xi)q 2 (x 2 ) 


All results are summarized in Theorem 1 in the text. 
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APPENDIX C 


NATURE OF COORDINATE SYSTEM CORRESPONDING TO K - 0} 4 0, k 2 4 0 

The following analysis will deal with the determination of the ex- 
plicit form of transformation (30). 

At the outset, the following observations are made. Suppose there 
are two coordinate systems (y-*-,y 2 ) and (xP-,x2) functionally related by 

y 1 = y 1 (x 1 ,x 2 ) 

If h^j and respectively, denote the metric tensors of the system 

(y^) and the system (x^), then - — - 

1J " 3/ SyJ ^ 

By differentiation of the above expression' coupled with other minor 
manipulations (ref. 20, p. 83), it is possible to show that 


ay 

dx 1 y 


: y 

ij) dx z 


rs J dx 1 dx<3 


where 


denote Christoff el- symbols based on g„„ and 


_ - 

hj_j, respectively. (See brief review at end of this appendix.) If the 

system (x^-,x2) is restricted to be a Cartesian system, all Christoffel 
symbols vanish and the above expression becomes 

ay 

dx 1 3x^ __ fiji dx l 


This system is then a system of second-order partial differential equa- 
tions for the transformation y i = y i (x 1 ,x 2 ). By reverting now to the 
notation employed in equation (30), the system can be written as the 
following set of first-order partial differential equations: 


3x7 11 


i.r 


i = 1,2 


l = 1,2 
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where Y 6tands for either or Y 2 . Evaluation, of the Christoffel 

symbols gives the following results: 


1 

11 

1 

12 

11 ) 

2 | 

12 ) 


1 d In (h*) 2 
2— 1X1 -° 

x d In (h*) 2 2 . 


"3xv 


£>(h?) 2 


2(h|)2 SX 2 

1 a m (h * 2 ) 2 

2 5x^ = 


a£X2 

b2 X f 


2 

/22 


1 d In (h |) 2 

2 " 0 


0 


2(hf) : 


S(b|)2 

3x7“ 


bffi 

a 2 X| 


The system of equations (Cl) can therefore be written as 


dY 

3x^ “ 


SY * 

" 2 

*1 

a 2 

X 2 - 

3x7 

“ b 2 

H 

^i 

^2 

?i 


v~ + 
*1 

X 2 


«l 

?2 


X 2 

X 1 

^2 

b 2 

X l* 

dx 2 “ 

' a 2 

X 2 ?1 


(C2) 
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It is possible to give particular solutions for and Y 2 which 

satisfy these equations and also satisfy the orthogonality requirement- 

dYl bYj_ Sy 2 SY 2 . - 

c% 3x^ + “ 0 — 


The functions 


where 


Y-L * cX x X 2 cos Id In X L - J In X ; 


Yo = cXtXo sin d In X n - ■= In X 


1 " d 


a 2 + b 2 


d 2 = 


represent a particular solution of the system (C2). The most general 
solution of the system can be therefore written as 

- n + °uV' " < C4 l 

where the are the elements of an orthogonal matrix and the T-\_ 

are constants. Equation (C4) merely represents a rotation and transla- 
tion of the coordinate system defined by (Y-^,Yg ) . 

In order to better ascertain the nature of the (X-^,Xg) coordinate 
system, let ^ 


p — cX-jXg 


0 = d In Xi - 


In X 2 


be a transformation of the (Xj_,X 2 ) system to a polar coordinate system j 
that is, in terms of p and 0, 

Y-j_ = P cos 0 

Y P = p sin 0 


Now consider the coordinate line X]_ = X£ = constant. The equation of 
this coordinate line is given parametrically by 

p = cX?Xo 

« i Cos) 

0 = dOnX^-ilnX 2 
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Equation (C6) can be expressed as 


e 9 = (xo) 4 ^ 1 /* 1 

(C7) 

Hence, from equations (C6) and (C7), 


p = c(X°) d V 0d 

(ca) 

or 


p = (const.) e -0d 



Similarly, the equation of the coordinate line X 2 - (const.) = is 
expressible as 

p = cCx®) 1 / 3 ^ 0 / 3 - (C9) 

or 

p = ( const.) e 0 / d 

The curves defined by equations (C8) and (C9) constitute a system of 
mutually orthogonal logarithmic spirals relative to the (p,0) system if 
the (p,0) system is interpreted as a system of polar coordinates in the 
plane. 


CHRISTOFEEL SYMBOLS 

A brief review of the definition of Christoffel symbols is presented 
here as an aid in following the analyses. 

Consider a surface in which an orthogonal coordinate system (x^-,x2) 
is embedded. In such a system the square of the differential of arc 
length is given by 

(ds) 2 = h 2 (dxr^) 2 + h|(dx^) 2 (CIO) 


hi =. h^x^x 2 ) 
h 2 = hgCx 1 ,x 2 ) 


when 
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The metric tensor g^j (i » 1,2) associated with (CIO) is defined by 

S 11 = h l J g 22 = ^2' g 12 = g 21 ~ 0 

The tensor g^ may be defined by 

gU , g 22 . ± g 12 , g 21 ,: 0 

Sll 822 

The Christoff el symbol ] [ (sometimes written r * “is then defined by 

(ij) Jk - 

j Z l , g ;k / 5g ik , 5g jk _ 5 Sjj\ 

*ij^ 2 \bx^ dx* Sx 15 - / 

Repeated indices imply summation. 
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APPENDIX D 


ANALYSIS: K = 0; fc, S 0, kg ^ 0. PROOF OF THEOREM 2 

By f oll owing the same procedure that was used in establishing lemma 
1, it can be shown that equation (24) becomes 


-1 


u 2N 


k 2 ( 1 + 2n + a-L a 9 J = 0 

U I 

As k 2 4 0 and a^ and ag are nonzero constants, it follows that 

U X = Cl U 2 (Dl) 

as in the previous case. 

From the definition of k^, 

1 ^ h l 

k n = r-4— = 0 


1 = ^ ^2 


Hence, 


f 5( x l) 


(D2) 


Recalling the definition of the Gaussian curvature of the surface 
and the assumption that K = 0 gives 



From equations (L2) and (D3), 


K 

(D3) 


d fl ^h 2 \ 

^ \h-j_ dx^y " 0 
*'• h^ ^ " h 2 k 2 = f e( x 2) 


or 


k 2 = 


f 6 (x 2 ) 


(D4) 
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At this point, it should he noted that kg cannot he constant for, if 
kg 3hg -- . 

it were, hp = ■ ? — 7 r and -5 — = 0, which, would imply that kp = 0. It 

W*2> d3! l - teJ 

is also follows that - — = 0, for if it is assumed that - — is zero. 


dx-|_ 

then, from equation (D4) , 




5kg 


f 6 (x 2 ) dh 2 

~zr~ 53T = 0 


/ >, ahp 

Now, either fg^Xg) = 0 or ^ — = 0 implies kg = 0, which violates 
the hypothesis. ^ 

Assume that — ^ ^ 0. Then, from conditions ( 3 ) and © and the 
Bxg 


proportionality of U]_ and Ug, 


g 2 = h 4 U x kg 


(D5) 


. a In g 2 _ d ^ U 1 ^ ^ k 2 

ax 2 - ax 2 ax2 


tf* 

© 

*s! 


Hence , 


s 111 k 2 a m g 2 5 m Up 


dx.? “ dxc 


axe 


(D6) 


Under the assumption, the right side of equation (D6) cannot he zero. 
Hence, the o.d.e. conditions arising from the combinations © and 0 
and © and © along with equation (Dl) yield 


a In kg 

-3Z~ = * 5 * 2*2 


(07) 


Similarly, it can he shown that 


a In k. 


•33cf " Vl 1 ^ 


(D8) 


Differentiating equation (D7) with respect to x^ and employing equa- 
tion (D4) yield 


a 2 In kp p* 

- b 5 ji. b 2 k 2 = 0 


(D9) 


npiri 
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in 

r- 

CD 


Now, differentiating equation (D8) with respect to x 2 and recalling 

that s 0 (as k-^ has "been assumed zero) give (with eq. (D9)) 

x 2 


S 2 In kc 


Skc 


5xg 


0 


&k‘ 


Thus, ^ 0 is incompatible witb kg ^ 0, k-^ = 0, and the only remain- 
2 

ing possibility is that kg is a nonconstant function of x^_ alone. 
Defining kg = fy(x]_) gives from equation (D4) 

f 6( x 2) 


hz 


(DIO) 


Finally, from the definition of kg, 

S In hg 

k 2 = h£ 3xj 




^ = 


S In hg 
^x^_ 


f 7(*l) 

fyCxq) 


(Dll) 


From these results. Theorem 2 follows directly. 

Because of the symmetrical nature of the equations , similar results 
can be obtained if it is first assumed that kg = 0 and kj_ ^ 0. Xn 
this case, 

kl ■ f s( x 2) 


hl- 


f 9( x l) 

f 8( x 2) 


and 


2 
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APPENDIX E 


PROOF OF THEOREM 4 

The four cases presented on page 29 will now be analyzed. 

(l) Case A: Ug = const. 

# With Uo constant, the following o.d.e. condition from (J)* and 


(D (set (45)) results: 


5 In hg/kg 

S In hg- 1 - 1 







. S 111 h 2~ Cll / k 2 n 

• • ^ =3 (J 

OX^ 



Hence,, 

— 


• 

hg" C11 = k 2 Pi(x 2 ) 

Pl( x 2) 

fo 

(El) 

In a similar manner from ( 4 )* and (8) 

* there result 


hj C±2 = k 2 <li(x 1 ) 

Ql( x l) 4 0 

(E2) 

Substituting according to equations (43) and (44) 

into equation 

(El) 


gives 


[ci 0 f i 3 (x 2 )e ^)] 2 - 011 = -i- re‘V- 2 ) 


3 - c nr f , ( )i 

. C 10 L r 15^ x 2 ; J 


2-cii 


( Cll -3)<P 

<p'e 


(E3) 


Taking the partial derivative of equation (E3) with respect to X]_ and 
assuming c^ ^ 3 give ^ — 
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Hence , 

(cn-3)cp 

= constant 

du 

. (cu-S)* 


or 

. ~ ( c i 2 ~5) 

9 = ln(au + b) 

(E4) 

A similar analysis 

using equations (E2), (42), and (44) yields 



, ~ ,-*-/( c i2~3) 

cp = ln(au + b) 

(E5) 


for Cj2 *£ ^r* 311 equations (E4) and (E5) it follows that 


C 11 ~ C 12 

If c -| i = 3, there results frcan equation (E3 ) that cp’ = constant and 
hence 

cp = au + b (E6) 

Next consider the case where one of the partial derivatives of U 2 
does not vanish and determine the corresponding form for cp. 

dU? , 

(2) Case B: y-=-sf 0; ^-=- so 

ox-l 7 ' dx 2 

The o.d.e. condition for and CD* gives 

d In Up 5 In hp 

- 5 ^ - c i3 -s- 5 w 

As U 2 is assumed to be a function of x^ alone, from equations 
(E7 ) and (43) 

r 2 ( x i) ■ =13 —5^— ’ “is? w 

= c 13 cp , (u)f ; [ 2 (x 1 ) 


(E8) 
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It follows from equation. (E8) that cp ' (u) => constant. Hence, 

cp = au + b 

An expression for U2 can he obtained from equations (E7) and (E8) by- 
solving 

a In U 2 , , , 

(const. ^ __ 

There results 

U 2 = (const . )e C ^-^^ - (E9) 

It can be shown that all remaining o.d.e. conditions are satisfied. 

dUo SUp 

(3) Case C: 5^ s °i ^ ¥ 0 

From symmetry considerations, this case is similar to case B, and 
the following results can he stated directly: 

9 =* au 4* h 

Ug = (const. )e C 15*13t x 2) 

BUo SUp 

(4) Case D: ^ 0; ^ 0 

The o.d.e. condition for (l) * and (5) * gives 

5 In Ug d In hg 

cix^ 

Hence, 

U 2 = h 2 16 q 3 (x2) ' (E10) 

Similarly, from the o.d.e. condition for ©* and (8)* 

Cl 7 

U 2 = h x p 3 (x 1 ) 


(Ell) 
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Equations (E10) and (Ell) give 


h2 16 = P3( x l) 

17 " 


(E12) 


Substituting according to equations (42) and 
assuming c-^g ^ c^y yield 


(43) into (E12) and 


( c i6 _c r7)<P 


[ g 10 f 15( x l)] 

tlO^S^)] 


•17 


1^7 = p^k^) C*u>) 


16 U 3 V 2 


From, equation (E13 ) , 


Since 


( c 16 “ © 3 . 7)9 = In P4(x x ) + In ^(xg) 



<p = au + b 


Hence , from equation (E13), 


P 3 (x 1 ) = (const. )e ClQf 12 Jf 12 (x 1 )J 17 


(EL4) 


(E15) 


and 


q 3 ( x 2 ) = (const. )e 


c 18 f 13 


(x P ) 




-C 


16 


where 


C 18 “ ^ C 16 _ C 17^ a 


c 19 “ ”^ c 16 " c i7) ac 2 
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Now consider the o.d.e. condition involving and 

5 In hg/Ughg d In hg 

= c 20 ^ 1 - 


d In h|" C20 /u 2 k £ 

■' 5 


= 0 


(E16) 


Substituting according to equation (E10) in (E16) gives 

a i4~ ( ° 2aKl6, A 2 


dxn 


= 0 


r -I d In hg d In kg 

•*•[2 - (c 20 + C 16 )J Sx£ ” 5 x£ = 0 

Therefore, from equations (43) and (44), 

[ 2 - ( c 20 + = 16 ^' + V ' = 0 

As cp' ^ 0, it follows that 


c 20 + c 16 = 3 — - 


From (4)* and a similar relation between constants results, and no 

new information is obtained. Remaining o.d.e. conditions are readily 
shown to be satisfied. 


Finally, consider the possibility c^g = Cj_j. Equation (E12) then 
becomes I 



and it follows at once that 

- C 2l[ f l2^1 ) ] ° 16 ‘ 

g 3 (x 2 ) » G 2lf f 13^ y -2^J 16 

At this point, the analysis for 9 follows the same pattern as given 
for case A, and the same equations for 9 are determined. Theorem 4 
follows at once. ' 
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APPEHDIX F 


ANALYSIS: Ui = Cpl^i ki f 0, kg a 0. PROOF OF THEOREM 6 

Four possible cases are distinguished as in previous analyses. 


Case A: Ug = constant 

Setting Ug = constant results in three o.d.e. conditions which 
must be satisfied. These conditions can be obtained from the list in 
set (45) and, after simplification, the terms that must be proportional 
are the following: 

* 1 9 In lu 

© 2k 2 " ET 


©* b 


^ ' 5x x 

-L S in k 2 / dh! 

note: s 0 as k]_ 


dx c 


■«) 


<D* k 2 


The o.d.e. condition for (D* and (§)* is 

*i S In k ? 

2k 2 - ^ " C 22*2 


Hence, either 


5 In kg 

5x1 


= 0 and Cgg = 2 


(FI) 


or 


Sk^ 1 


C c 22 " 2 ) h l = ^ c 22 ’f 2 ) 


Also from and 


there is the possibility that either 
5 In kg 


5xc 


= 0 


(F2) 


(F3) 


or 


h 2 = c 23 


5k 

5xc 


-1 


(F4) 
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If equations (Fl) and (F3) hold, then kg = constant and 

h 2 = pJCx^q^Xg) 

Pl( x l) 


If it is assumed that equation (F4) is valid, there is obtained upon 
differentiation 

dh.g S^k-l. 

^ = C 23 

However, from equation (Fl) or (F2) it can be shown -that s 0. 

Sh 2 ^ 1 

But this would 'imply that ^ = h^hgkg =0, which violates the hypo- 
thesis. Hence, the one remaining possibility is 

S In ki S In ko , 

—xr-*- ¥° 


Therefore, 


k 2 = PiCx-L) 


From equation (F2) it then follows that 

_ 1 Pl( x l) 

1 2 “ C 22 p2( Xi ) 


3 In he 


Pl( x p) 


*x-L 2 - c 22 P*(x i ; 


r -iV( 2-c 22 ) 

^2 = [Pl( x l)J <ll( x 2) 
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dU 2 dU 2 

0assB: S£#°>S[ s0 

By an. analysis similar to that given for case A, it can he shown 
that h-j_j h 2 , kg, and U 2 mus '* : have the following forms: 


U 2 = (const. 

Pl(xi) 


hi = 


Pl( x l) 


q.i(x 2 ) 
2 " p^(^T 

= Pi( x i) 


du? , dU 2 

Case C: — =k 0} •*-— =0 

OXi OXg 


The solutions for hi, hg, kg, and U 2 are U 2 = ^Pi( x i) + const .J 

— h 2 

with ^ and h^ defined hy equation (P5) or ^ = (const.) — y 

with k 2 , hi, and h 2 defined hy equations (P7), (P8), and (F9), 
respectively. 


dUc 


£Uc 


Case D : ■ zjk 0 and ^ 0 

The solutions for h l> hgj kgj and Ug are 

u 2 = (const. Jp^Cx^e *^ qi ^ 2 ^ 2 (n ^ 0) 


h-, = - 


pjC x i) 

p 2 ( x l) 

. _ 

^ PltXp) 

kg = Pi(xi) 
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As mentioned previously, the case where kg = 0 gives similar re- 
sults, with the roles of and x 2 interchanged (from symmetry con- 

siderations). Theorem 6 readily follows. 
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